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Abstract 

We show that flat spectrum of small perturbations of field(s) is generated in a 
simple way in a theory of multi-component scalar field provided this theory is confor- 
mally invariant, it has some global symmetry and the quartic potential is negative. 
We suggest a mechanism of converting these field perturbations into adiabatic scalar 
perturbations with flat spectrum. 

1 Introduction 

An elegant mechanism of the generation of (almost) flat spectrum of scalar perturbations is 
provided by inflation [T] . Yet it is legitimate to ask whether the spectrum of the Harrison- 
Zeldovich type can emerge in a different way. This question is of interest, in particular, 
from the viewpoint of alternatives to inflation — pre-Big Bang scenario p] (for a review see 
Ref. [3]), ekpyrotic/cycling models [H [5] (for a review see Ref. [6]), starting-the-Universe 
picture [7], etc. One possibility is to make use of scalar fields with negative exponential 
potentials [H [9l [TOl [11] (see also Ref. [12] )• In this way one indeed obtains almost flat 
spectrum of scalar perturbations in ekpyrotic models, with overall homogeneity and isotropy 
of the Universe preserved during the contracting phase because of the stiff effective equation 
of state [13]. In concrete models of this type the contracting solutions are unstable [HI [15], 
so that one has to fine tune initial data in our part of the Universe. It has been argued, 
however, that such a fine tuning may be automatic, in a sense, in a cyclic scenario with long 
periods of dark energy domination [16j . 

In this paper we suggest another mechanism for generating the Harrison-Zeldovich spec- 
trum. The idea is to relate scale invariance of this spectrum to conformal invariance of field 
theory. We will see in Section [2] that this can be done by introducing conformal scalar field 
as a source of entropy perturbations, which are converted into adiabatic perturbations at 
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later epoch. The three basic ingredients are that (i) the theory possesses some global sym- 
metry, (ii) the quartic scalar potential is negative and (iii) there is long enough stage of the 
cosmological evolution prior to the conventional hot Big Bang epoch. At that early stage, 
the classical scalar field rolls down its potential along the late time attractor. In this regime, 
the vacuum perturbations of field(s) orthogonal to the radial direction — the direction along 
which the classical field rolls down — get amplified and freeze out having the spectrum which 
is automatically fiat. Hence, getting fiat spectrum of the field perturbations is reasonably 
simple. 

Implementing this mechanism to obtain adiabatic matter perturbations of correct ampli- 
tude is somewhat less straightforward. We propose the corresponding scenario in Section [3J 
At the end of the rolling stage, both conformal invariance and global symmetry are assumed 
to be broken, and the scalar potential is assumed to have a minimum or a set of minima, 
see Fig. [1] If the scale of global symmetry breaking is small enough, the situation becomes 
analogous to the axion misalignment picture (for a review see Ref. [IZj), the analog of axion 
being the orthogonal, pseudo-Goldstone field(s) corrseponding to broken global symmetry. 
Unlike in the case of axion, the interactions of the pseudo-Goldsone field with matter are 
not negligible, so this field plays the same role as curvaton in infiationary models. It gets 
frozen at a slope of the potential, later on it starts to oscillate and eventually these oscilla- 
tions convert into conventional particles. In this way the adiabatic perturbations with fiat 
spectrum are produced. 

In Section H] we discuss a subtle point of our model. Namely, at the stage when the 
radial field rolls down, this field develops its own perturbations. Their growth is fast, and 
at some point the linearized theory breaks down. We argue, however, that the results of the 
linearized theory concerning the perturbations of the pseudo-Goldstone field remain valid, 
and the fiat spectrum of adiabatic perturbations is generated indeed. 

We end up in Section [5] with concluding remarks. 

2 Getting flat spectrum 

To begin with, let us consider massless scalar field conformally coupled to gravity and 
evolving in spatially fiat 4-dimensional FRW Universe. Let us assume for the time being 
that there is exact global symmetry, which we take for definiteneness to be U{1) (though 
the argument is valid for any compact symmetry group). Hence, the field is comlex, and 
the global symmetry is — > e*"0. Conformal invariance allows for quartic scalar potential, 
which we assume to be negative, 

V{<P) = -h'\<p\' 
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Figure 1: The scalar potential is negative quartic and respects global symmetry (f/(l) here) 
at relatively small 0, and has one or more minima at large 0. Dashed line shows the evolution 
of the field (dot with arrow) at the rolling stage. At that stage the perturbations (shown by 
double arrow) of the pseudo-Goldstone field are developed. At the hot Big Bang epoch, the 
pseudo-Goldstone field transfers its energy to hot matter, and its perturbations are converted 
into adiabatic perturbations. 

(modulo an additive constant). As usual, one defines the field x by 

^ ^ X ^ Xi + iX2 
a a 

Then the action for x in conformal coordinates is 

S^ = j d'xdv [d.x^d'^x + h'lxl'] , 

where indices are raised by Minkowski metric. 

Let us consider spatially homogeneous background field. The field equation is 

~x" + 2h'\xW = 0, (1) 

where prime denotes d/drj. In terms of the radius and phase, x = P^'^, one of the equations 
is the conservation of current. 
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Hence, as the value of p increases, the phase 6 freezes out, and the evolution proceeds along 
the radial direction. Without loss of generality we take x to be real at that time, x = Xi- 
Then the first integral of the equation of motion is x'l^ ~ h'^Xi = ^ = const. Hence, the 
solutions approach the asymptotics that is independent of e, 

X-.W = . (2) 

where r/* is a constant of integration. 

Now, let us study field perturbations in the orthogonal direction, 8x2- They obey the 
linearized equation, in momentum representation, 

{5x2)" + k''5x2 - 2h\l5x2 = . (3) 

At early times, when A;(?7* — r^) ^ 1, the second term dominates, and 8x2 oscillates like free 
scalar field in Minkowski space-time. At later times the third term dominates instead. So, 
the perturbation "exits the horizon" , but now the "horizon" is due to evolving xi (hereafter 
we use quotation marks to distinguish this "horizon" from the true cosmological horizon). 
Outside the "horizon", i.e., at /c(?7* — 77) -C 1, the perturbation 8x2 evolves in the same way 
as the background xi- This is clear from the fact that at A; = 0, eq. ([3]) is the linearized 
equation ([1]). Now, alongside with the solution sq. ([I]) has a solution e^^Xiiv)^ where a 

is a real constant. For small a the latter solution is (xi + i(^Xi)i ^i-nd the second part is small 
perturbation, which is precisely 8x2- So, if the perturbation 8x2 oscillates with amplitude 
c{k) at early times, it behaves at late times as 

8x2 = Cc{k)xM = C'-^^^, (4) 

where the factor k~^ is evident on dimensional grounds, C and C are time- independent, 
and C is independent of k. This is precisely the enhancement of perturbations needed for 
getting fiat spectrum of 8x2- We stress that this property is entirely due to the symmetries 
of the model: the behaviour ([2]) of the background is a consequence of conformal invariance, 
whereas the behaviour (jlj) of perturbations is a consequence of the global U{1) symmetry. 
In more detail, eq. ([3]) reads 

(8x2)" + kHx2 - , ^ ,, ^X2 = . 

This is formally the same equation as the equation for perturbations of minimally coupled 
massless scalar field in de Sitter space-time. Hence, the spectrum of 8x2 is the same fiat 
spectrum. We are interested in the solution that behaves at early times as 

' (27r)3/2^ 
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Modulo overall constant phase, this solution is expressed through the Hankel function, 



At k{rj^ — Tj) <^ 1 this solution is 

in accord with (jlj). Assuming that the quantum field 5x2 is originally in the vacuum state, 
we obtain its power spectrum outside the "horizon" in the standard way. The spectrum is 
fiat, and its amplitude is 



As promised, generating fiat spectrum of the field perturbations is fairly straightforward. 

We end up this section by noting that the mechanism requires long time during which 
the classical field rolls down its potential. Indeed, let rji be the time at which the field 
starts to roll down, and rjf the time at which rolling down terminates. Then the existence 
of the regimes ([5]) and ([6]) requires that 

These inequalities should hold for all scales of interest, say, for present momenta (we set the 
present value of the scale factor equal to 1) 

kmin ~ (1 Gpc)-^ < A; < krnax ~ (100 kpc)-^ . (8) 

Therefore, the duration of the rolling stage should be large. 

This excludes the possibility that the mechanism works at the hot stage of the cosmological 
evolution: there is simply not enough time. On the other hand, there may well be enough 
time in the ekpyrotic or starting-the-Universe scenaric0. 



^Note that the inequahty ^ means that the comoving size of hght cone originating from t] — rji and 
measured at rj — rif exceeds the comoving size of the visible Universe. Hence, the mechanism can work only 
in cosmological models which solve, at least formally, the horizon problem. 
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3 Implementing the mechanism 



In this Section we suggest a way to implement the mechanism and generate adiabatic per- 
turbations from the perturbations 502- Let us assume that as the classical field 0i = Xi/*^ 
approaches the region 0i ~ /, conformal symmetry gets broken and the potential along the 
radial direction has a minimum at 0i = /, where / is some high energy scale. We assume 
that 0(2) symmetry is broken as well, but that the scale of that breaking is small compared 
to /, so we neglect the effects of 0(2) symmetry breaking for the time being. Then the 
classical field 0i oscillates near / and eventually settles down to 0i = / (say, by producing 
conventional particles). We require that this process occurs at the time when the energy 
density of the field 0i is much smaller than the total energy density in the Universe. The 
reason for this requirement is that the field 0i develops its own perturbations. In terms of 
perturbations 5xii the linearized equation for this mode is 

{5x1)" + k-'dxi - !^)2 '^Xi = . 

The solution which is negative frequency at early times is again expressed through the Hankel 
function, but with another index. 



At k{ri^: — rj) <^ 1 this solution is 

v(-) = I I (10) 

7rr(-3/2)A;5/2(ry^_ry)2 • ^ > 

This behaviour again has simple interpretation: if xi is a classical solution to the field 
equation ([T]), then 

^Xi x'l (11) 

is a solution to the linearized classical field equation for A; — > 0. Hence 6xi oc (r^* — 1])^ on 
super-" horizon" scales; the delpendence on k again follows from dimensional argument. In 
any case, eq. (ITO!) implies that the spectrum of perturbations of the field Xi is red. 

These perturbations should not leave any trace in the late Universe; the simplest way to 
achieve that is to impose the above requirement. 

The perturbations of the field 0i may become so strong that they ruin the linearized 
analysis. To be on the safe side, we impose here the requirement that the field 0i is approx- 
imately homogeneous over the whole visible Universe at the time when the scales of interest 
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exit the "horizon" (see, however, Section |4]). This is sufficient for the hnearized theory to 
be vahd in our patch of the Universe at that time; perturbations S(f)i of longer wavelengths 
correspond to unobservable shift of the parameter rj^ in the whole visible Universe, see (fTTj) . 
Our requirement gives 

where is the co moving size of the visible Universe, cf. ([8]), and i^max is the time of the 
"horizon" exit of the shortest interesting modes. The value of rimax is determined by 

kmax ij]* '^max ) ~ 1 • 

Making use of ([2]) and (fT2|) we obtain a constraint on the coupling constant 

h^^^ 10-4 . (14) 

We will further discuss the non-linearity issue in Section HI 

Let us now turn to the field 02- In fact, after the field 0i settles down to /, it is more 
appropriate to speak about the phase 9 of the original field 0. We keep the notation 02, with 
understanding that = 02//- 

Let us ffist estimate the amplitude of perturbations 502 at the time rjf when the back- 
ground field 01 reaches /. At that time one has 

XiM~^^(;^-a/-/, (15) 

where the subscript / refers to the time rif. From ([7]) we obtain 

We now have to convert these entropy perturbations into adiabatic ones. 

Let us recall that the global 0(2) symmetry is broken near |0| = /. We have in mind 
axion-like picture, in which the energy scale of this breaking is smaller than /. Importantly, 
the phase field 02, as any (pseudo-)Goldstone field, minimally couples to gravity fT8]. From 
this point on, the discussion is parallel to that in the curvaton model. 

Generically, 02 = is not a minimum of the scalar potential, so the radial evolution ends 
up at a slope of the potential for the pseudo-Goldstone field 02. Let us assume that at the 
time rjf, the Hubble parameter exceeds the mass parameter of this pseudo-Goldstone field. 
In this case the field 02 stays close to zero for some tim^, then rolls down to the nearest 
minimum of its potential, oscillates there and eventually the oscillations decay into usual 



^This is true both for expanding Universe and for contracting Universe dominated by matter with stiff 
equation of state, w > 1. 
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particles. We assume that the decay happens at the hot Big Bang stage of the cosmological 
expansion, when the Universe is either radiation-dominated, or temporarily dominated by 
the oscillating field 02 itself. Let M denote the distance to the nearest minimum of the 
potential for the filed 02, so that the number of minima of the scalar potential is ~ f /M- 
Then at the beginning of oscillations and until their decay, the perturbations in the energy 
density of the pseudo-Goldstone field are 

This leads to adiabatic perturbations with fiat spectrum and the amplitude 

A~r^^^~r4-, (16) 
ajM 2'kM ' ^ ' 

where r is the ratio of the energy density of the field 02 to the total energy density at the 
time the oscillations decay. 

The pseudo-Goldstone picture with a single minimum of the potential for the phase would 
correspond to M ~ /. If we insist on the constraint (fT^ . the amplitude (fT6!) of the resulting 
adiabatic perturbations would be too small in that case (barring fine tuning). On the other 
hand, with ~ f /M 3> 1 the required amplitude A ~ 10~^ can obtained without fine 
tuning of other parameter^. As we discuss in the end of Section HJ the constraint f|T3l) . and 
hence f|T^ may in fact be unnecessary, so the mechanism may work for A^ = 1 as well. 



4 Generalizing to non-linear radial field perturbations 

Let us come back to the non-linearity issue. Since the perturbations (fT2l) increase in time 
faster than the background field ([2]), the evolution in the radial direction may become non- 
linear at some stage. Naively, one would think that the analysis of Section [2] goes through 
only if this does not happen, i.e., if the perturbations 50i remain small up to the time rjf. 
Such a requirement would push the parameters of the model to a very contrived range: the 
coupling contant h would have to be extremely small, the field value / would have to be 
super-Planckian, and the number of minima of the scalar potential, A^ ~ f 1^-, would have 
to be very large. This is discussed in Appendix. However, we argue that non-linearity of the 
radial field perturbations does not spoil the result (fT6|) . 

•^This discussion assumes that after the rolhng stage, the value of the phase of the field (f) is generic. 
Alternatively, one could make use of the anthropic argument for explaining why the field 02 ends up close to 
the minimum of its potential, i.e., why M ^ /, so that the amplitude of the matter perturbations is large 
enough, cf. Refs. [HI [20]. With this line of reasoning, having large number of minima of the potential for 
the pseudo-Goldstone field is certainly unnecessary. 
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The argument is as follows. Once the inequality (fT3|) is satisfied, all interesting scales 
exit the "horizon" when the linearized theory is still valid. Let rjun be some moment of time 
when the linearized theory is valid, but all interesting scales are already super-" horizon" . 
At that time, 502 niay be considered as classical random field whose spatial gradients are 
negigible. According to ([7]) and ([2]) this field is small, 502/01 ~ h, and it evolves in the same 
way as 0i, 

^Mv) = —1 — ■ (17) 

The field 0i(x, 77) = 01(77) + 50i(x, 77) is also classical random field which is homogeneous 
inside the "horizon" . At late times this field becomes inhomogeneous, but on super-horizon 
scales only. The point is that the evolution of the entire system proceeds according to the 
homogeneous classical field equation. Therefore, the symmetry argument presented before 
eq. (jl]) tells that 502 is still given by f|T7|) . but now with 0i = 0i(x, 77). Note that 502 remains 
small. The field 0i reaches the value / in the end, and at that time 

(pliVlin) 

This coincides with the result of the linearized theory. 

Another way to phrase this argument is to write the rolling field 0i, with long-ranged 
perturbations included, as 

0i(x,r/) 



/7,a(7/)[7/,e//(x) -T]] ' 



Comparing this expression with fll2p at the time when the linearized theory is still valid, we 
find that 57/* (x) = rj^ eff{^, v) ~ V* is long ranged random classical field whose fiuctuation is 

3h 

A, 



2nk 



At late times, when 57/* > (7/* — 7/) , the dependence of 0i (x, 7/) on 57/* becomes non-linear, 
but the expression f[T5]) remains valid. The point is, again, that the expression f[T7l) with 
01 = 01 (x, rj) is a solution to the linearized field equation for 502 at all times. To see this 
explicitly, let us show that |V502| <^ 502, so that spatial variation of 7/*e// can be neglected 
in that field equation, and the argument before eq. (jlj) indeed applies. Using (fTTIl we write 

= |V«,.| . (19) 



Now, V577* is random field with fiuctuation 
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so the ratio (fT9|) is indeed small. Since the final state of the super-" horizon" perturbation 
502 is independent of the parameter 77*, we see that non-linearity of the evolution of 0i, 
encoded in the non-linear dependence of 0i on ^r^^,, is irrelevant for the spectrum of 5(j)2- 

The latter way of presenting the argument suggests that the condition (IT^ is in fact 
unnecessary. The only requirement that remains is that /i -C 1. This ensures that the 
perturbations 5(f)i become non-linear long after they exit the "horizon", and the spatial 
variation of ?7*e// is irrelevant in the field equation for 502- Hence, the parameters of our 
model are constrained very weakly; in particular, there may exist just one minimum of the 
potential for the pseudo-Goldstone field, i.e., M ~ /. 

5 Conclusions 

We conclude this paper with a few remarks. First, the mechanism for converting the per- 
turbations of 02 into adiabatic perturbations, suggested in Section [2l is most probably not 
unique. For example, one can imagine adding other fields which interact with the pseudo- 
Goldstone field 02 so that the energy density of these fields, after 0i relaxes to /, is linear 
in 502. In that case the simple formula flTB]) would not be valid. 

Second, the resulting spectrum needs not be exactly fiat. One expects that for conformal 
symmetry slightly broken during the stage of rolling 0i, the spectrum becomes tilted. The 
tilt may be present also if the global symmetry {U{\) in our example) is explicitly broken 
at that stage. How strong is the variation of the tilt with fc, and whether there are other 
properties of the spectrum in these cases, remains to be understood. 

Third, if the coupling constant h is not very small, one would expect sizeable non- 
Gaussianity of the perturbations. Also, there may be an interplay between the perturbations 
501 and 502 which may result in peculiar properties of the adiabatic perturbations. It would 
be interesting to see whether or not the latter may help discriminate our mechanism from 
others (say, from infiationary ones). 

Fourth, as in many other cases, the generation of scalar perturbations in our model is 
unrelated to the generation of tensor modes. Thus, there is no reason to expect sizeable 
gravitational wave background in our scenario. 

Finally, we left aside the issue of the initial conditions for the field 0. We simply assumed 
that this field is spatially homogeneous and starts rolling from near the top of its quartic 
potential. Whether these initial conditions may emerge naturally in the context of some 
cosmological scenario remains unclear. 

We hope to address these and other issues in future. 

The author is indebted to D. Gorbunov, A. Khmelnitsky, E. Nugaev, G. Rubtsov and 
I. Tkachev for helpful discussions. This work has been supported in part by Russian Foun- 
dation for Basic Research grant 08-02-00473. 



10 



Appendix 



In this Appendix we discuss what kind of constraints on the parameters one would obtain 
by insisting that the entire evoluton of both perturbations 6(f)i and 502 occurs in the hnear 
regime, so that the analysis of Section [2] is literally valid. The strongest constraint comes 
from the requirement that the perturbations of the radial field 6(j)i are small by the end of 
the regime of rolling 0i, 

6MVf)-^xAVf)/(^f'^f- (20) 

Let us pretend that this inequality should be valid for length scales up to the present hori- 
zon size /Cmin) perturbations of longer wavelengths correspond to unobservable shift of the 
parameter rj^, in the whole visible Universe. 

Making use of (fT5|) and (fT2l) one finds that the condition (1201) gives 

^ » h'f . (21) 
a/ 

Given the very small value of kmin of the order of the present Hubble parameter, the latter 
inequality impies that the coupling constant h is extremely small. Indeed, assuming that 
the Universe did not expand significantly before the beginning of the conventional hot Big 
Bang stage (no inflation) and ignoring weak dependence on the effective number of degrees of 
freedom g^,, one finds that a/ > To/Tmax, where Tmax and Tq are the maximum temperature 
in the Universe and the present temperature, respectively. Hence, the constraint (12T]) reads 



h'f « ^T^ax - lO-''T^ax < . (22) 

The constraint (pTIl implies also that 

f:^Mpi. (23) 

This comes out from the estimate of the temperature Tend at which the oscillations of the 
pseudo-Goldstone field 02 are converted into usual particles. At this temperature the per- 
turbations 502 are transformed into adiabatic ones, and one requires that T^nd ^ 100 GeV in 
order that the baryon asymmetry and dark matter may be thermally produced afterwards. 
To estimate Tend in our model, we notice that the oscillations of 02 begin at the time tosc 
when H{tosc) ~ where /x^ is the curvature of the potential for the pseudo-Golstone field. 
Again ignoring the dependence on g^, we obtain the energy densities of hot matter and the 
pseudo-Goldstone field at that time, 

Prad{tosc) ~ T'^{tosc) ~ , (^osc) ~ /i^M^ . 
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The oscillations end up when Prad/ P<p2 = ^ < 1, so one finds 

Tend (^(tosc) 1 ^ P(h2^osc) 

rsj ^ 

Tqsc 0,{tend) ^ Pradij^osc) 



Hence 



Making use of (fT6l) we find 



Tend-^:^^pMpi 



Finally, using (!22|) and inserting A ~ 10 ^ we obtain 



Tew « ^ GeV . 



Since p"^ is the curvature of the scalar potential, one has p < Mpi. Therefore, high enough 
value Tend ^ 100 GeV can be obtained only for super-Planckian values of /. Together with 
( 122|) this implies that the coupling constant is extremely small. In turn, the smallness of h 
and eq. flTBl) mean that the correct amplitude of perturbations can be obtained only for huge 
number of minima of the scalar potential, N = M/ F. 
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